Dissipative Diamagnetism — A Case Study for Equilibrium and Nonequilibrium 
Statistical Mechanics of Mesoscopic Systems. 
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Using the path integral approach to equilibrium statistical physics the effect of dissipation on 
Landau diamagnetism is calculated. The calculation clarifies the essential role of the boundary 
of the container in which the electrons move. Further, the derived result for diamagnetization 
also matches with the expression obtained from a time-dependent quantum Langevin equation in 
the asymptotic limit, provided a certain order is maintained in taking limits. This identification 
then unifies equilibrium and nonequilibrium statistical physics for a phenomenon like diamagnetism, 
which is inherently quantum and strongly dependent on boundary effects. 
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An unconventional approach to statistical physics, 
which may be referred to as the Einstein approach, in- 
volves the derivation of equilibrium results from the long- 
time limit of time-dependent equations 0. Specifically, 
a set of Langevin equations (or their equivalent in the 
phase space, called the Fokker- Planck equation), with 
built-in detailed balance conditions, can naturally yield 
asymptotic results that can be independently calculated 
from the Gibbs ensemble idea of statistical physics. The 
underlying concept is physically appealing because not 
only does it sidetrack the issue of ergodicity, which is 
assumed at the outset in the Gibbs prescription, it also 
connects directly to experimental measurements, which 
necessarily involve time-averages. In this centenary of 
Einstein's annus mirabilis it is momentously appropriate 
to assess the validity and usefulness of this approach to 
statistical physics, that relies on the central paradigm of 
Brownian motion 

Given this motivation we want to further explore the Ein- 
stein approach in this Letter by going beyond the classical 
into the quantum domain. The phenomenon of interest 
happens to be intrinsically and essentially quantum me- 
chanical — it relates to the issue of diamagnetism exhib- 
ited by a collection of electrons subjected to an applied 
magnetic field. Diamagnetism is an enigmatic subject 
in that not only does it require a quantum treatment, 
as provided by the landmark work of Landau , but it 
also needs a careful analysis of the boundary of the con- 
tainer in which the electrons are constrained to move. As 
has been discussed lucidly by Van Vleck Q, the bound- 
ary electrons exactly cancel the contribution of the bulk 
electrons, in classical physics, leading to the celebrated 
Bohr- Van Leeuwen theorem [5|. However this cancella- 
tion is incomplete in the quantum regime, because as 
Peierls points out it is the boundary electrons which 
have the "skipping orbits" that yield the edge currents, 
familiar also in quantum Hall effect Q, which make an 
essential contribution to diamagnetism. A few years ago, 
we have examined the question of Landau diamagnetism 
in a dissipative and confined system ||. 
The following issues were addressed in I: (a) the approach 



to equilibrium of a quantum dissipative system, the anal- 
ysis of which brings out the subtle role of boundary elec- 
trons, (b) the effect of dissipation on Landau diamag- 
netism, an equilibrium property, (c) quantum - classical 
crossover as the system transits from the Landau to the 
Bohr- Van Leeuwen regime as a function of damping and 
(d) the combined effect of dissipation and confinement on 
Landau diamagnetism, the latter arising from coherent 
cyclotron motion of the electrons. The item (d) is par- 
ticularly relevant in the context of intrinsic decohcrence 
in mesoscopic structures in view of heat bath induced 
influence 0, 0, [^}. Dissipation was incorporated in I 
with the aid of a quantum Langevin equation, driven by 
a systematic Lorentz force, that can be derived from an 
underlying Hamiltonian in a system-plus-bath formula- 
tion in which the bath degrees of freedom are integrated 
out In the infinite past the bath is assumed to be 

in thermal equilibrium such that the fluctuations of its 
degrees of freedom are governed by quantum statistics. 
Thus, detailed balance conditions are automatically ex- 
pressed through a 'fluctuation-dissipation' relation that 
relates the noise spectrum to the damping term in the 
quantum Langevin equation. 

The starting point of I as indeed in this Letter is the 
Feynman- Vernon [T^ | Hamiltonian for a charged particle 
e in a magnetic field B: 
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where the first term is the Darwin [13j term represent- 
ing a confining potential, p and x are the momentum 
and position operators of the particle, Pj and Xj are the 

corresponding variables for the bath particles, and A is 
the vector potential. The bilinear coupling between x 
and Xj as envisaged in Eq. (1) has been the hall mark 
of the Caldeira-Leggett approach to dissipative quantum 
mechanics 113 • Assuming the B field to be along the 
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z-axis, all the vectors in Eq. (1) can be taken to lie in the 
xy-planc. From the quantum Langevin equation , derived 
from Eq. (1) by following the steps mentioned above, the 
nonequilibrium or time-dependent magnetization along 
the z-axis, M z (t) is computed in I. It is important to 
note that the Landau answer for the magnetization, in 
equilibrium, ensues from M z (t) only by following the lim- 
iting procedures in a specific order, viz; by first taking t 
— > oo and then setting ujq — > 0. If these two limits are 
interchanged one ends up with a piece of the Landau an- 
swer that misses out the boundary contribution. 
Having laid down the background to the myriad perplex- 
ing issues concerning diamagnetism we pose and answer 
the following question in this Letter. Should we not be 
able to calculate the equilibrium magnetization directly 
from Eq. (1) by following the usual Gibbsian statistical 
mechanics in which all the terms in Eq. (1) are treated 
on the same footing and there is no separation between 
what is a system and what is a bath? If the answer to 
this question is in the affirmative and the resultant mag- 
netization matches with the result derived in I in the 
'equilibrium limit' that would indeed lend the Einstein 
approach yet another foundational basis. 
The energy eigenvalues for the Hamiltonian in Eq. (1) 
have been computed by Hong and Wheatley 0|. How- 
ever our method of calculation is based on the functional 
integral approach to statistical mechanics which we find 
to be the most convenient tool for stud y in g charged par- 
ticle dynamics in a magnetic field [l7 l 111 H I |20L l2lj . 
The Euclidean action corresponding to the Hamiltonian 
in Eq. (1) can be written as : 
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where the subscripts S, B and I stand for 'system', 'bath' 
and 'interaction' respectively. The corresponding La- 
grangians are enumerated as: 
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x(r) 2 + w 2 f(r) 2 - oj c (x(t) x x{t)) z , (3) 



where uj c = j§-, is the cyclotron frequency, 
N 1 



(4) 



where the Bosonic Matsubara frequencies v n arc given by 
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The system-part of the action in terms of Fourier com- 
ponents is: 

A " = y ^E[(^ + w o)(^K)-^K)) 

n 

+uj c u n {x{v n ) x x*(v n )) . (9) 
In deriving Eq. (9) we have used the identity: 
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Following the detailed treatment given by Weiss [2Cj , the 
combined contributions of the bath and the interaction 
terms to the action can be written as: 
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Introducing the spectral density for bath excitations as: 

N 

J ( w ) = o H m ^ 3( K w ~ (I 3 ) 



we may rewrite 
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Now combining Eq. (11) with Eq. (9), the full action 
can be expressed as: 

A e = — Upy^^U 2 + U)l + Vnj(v n ))(x(v n ) ■ X* (l> n )) 
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+uj c i> n {x(v n ) x x*(v n ))\, 
where the 'memory-friction' is given by 

J (uj) v n 



(15) 



(16) 



We introduce now imaginary time Fourier series expan- 
sion of system variables and bath variables as follows: 



Note that x{v n ) is a two-dimensional vector 
(x{y n ),y(v n )). Introducing then normal modes: 
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Eq. (15) can be rewritten in a 'separable' form: 



{z+(v n )z* + {v n )) 

+(i>l + ul + v n *f(v n ) - iu c v n ) 
(*_("»)*-("»))]• (18) 
The partition function is then given by: 
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In view of Eqs. (8) and (16) the Helmholtz Free energy 
T can be deduced from Eq. (19) as 

(3 V 2tt 
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where the first term is independent of the magnetic field 
and owes its existence purely due to the Darwinian con- 
straining potential. Equation (20) contains all the ther- 
modynamic properties, the most important of which is 
the magnetization given by the negative derivative of T 
with respect to B : 
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Equation (21) identically matches with the asymptotic 
[t — ► 00 ) limit of the expression obtained by Li etal |22l | 
from a quantum Langevin equation formulation. Further, 
in the so-called ohmic dissipation model for which 

J(to) = M-yu, (22) 

the expression (21), upon using the identity: 
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also yields the asymptotic result of I, for oj = (cf. Eq. 
(19) of I). The ohmic case is relevant for electron-hole 



excitations in a Fermionic bath whereas the non-ohmic 
case applies to a phononic heatbath [2(j ■ 
Equation (21) embodies several tantalizing results which 
deserve special comments: (1) The diamagnetization is 
one of the rare equilibrium properties which depends di- 
rectly on the damping parameter 7. Seldom is dissipa- 
tion discussed in text books within the realm of what 
we call equilibrium statistical mechanics, based on the 
Gibbs ensemble. The fact that 7 is a measure of dissipa- 
tion has been amply demonstrated in I, wherein we had 
shown how by increasing 7, M changes from the Landau 
to the Bohr- Van Leeuwen expressions — an example of 
coherence-to-decoherence transition in an open quantum 
system j^. (2) Diamagnetism as a material property is 
seen to be situated at the crossroads of thermodynamics 
and transport phenomena. The thermodynamic nature 
of the property is rooted on its being able to be calcu- 
lated from the free energy, as shown here. On the other 
hand, diamagnetism, like the Drude conductivity p3 |. is 
also based on transport mechanism in that it is related 
to the expectation value of the operator (r x v) (see I) . 
Because the velocity v appears explicitly, dissipative dia- 
magnetism naturally connects to the fundamental fric- 
tional material property, viz. resistance, in view of the 
fact that 7 _1 is related to the Drude relaxation time |25j . 
Again we are not aware of any other phenomenon that lies 
at the juxtaposition of thermodynamics, which is derived 
from a partition function and transport, that is usually 
treated in kinetic theory. (3) Normally, in statistical me- 
chanics, a thermodynamic limit is taken as a result of 
which surface contributions to bulk become irrelevant. 
However, for diamagnetism the surface enters crucially, 
as argued above; even though, there are fewer surface 
electrons than in the bulk, their contribution to the op- 
erator r in (r x v) is substantial. A remarkable feature 
of diamagnetism is the need to first calculate the magne- 
tization in the thermodynamic limit and then switch the 
boundary off i.e. by setting loq = 0. One related issue is 
the environment induced dissipation which happens to be 
a ubiquitous attribute of a mesoscopic system. Addition- 
ally, because for a mesoscopic system surface effects are 
non-negligible, the present study has a bearing on our un- 
derstanding of mesoscopic structures. While points (1), 
(2) and (3) connote to thermal equilibrium we want to 
now make a few remarks on the significance of our results 
for the approach-to-equilibrium, in the present context: 
(4) usually this question is discussed in a system-plus- 
bath approach, within a master equation for the density 
operator. The subject of quantum optics is replete with 
such approaches wherein the interaction between the sys- 
tem and the bath is assumed weak and is consequently 
treated in the socalled Born-Markov approximation [26| • 
The result is, although the approach to equilibrium does 
depend on relaxation parameters such as damping the 
equilibrium results themselves are independent of such 
parameters. Thus the density operator approaches a 
Boltzmann distribution characterized by the Hamilto- 
nian for the system alone. In contrast, the presently 
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derived dissipativc diamagnetization, which can also be 
computed from the nonequilibrium method of I, does de- 
pend explicitly on damping, as has been also emphasized 
under point (1) above. The reason is, like in the much 
studied problem of quantum dissipation of a harmonic 
oscillator [23] , the system- bath coupling is so strong that 
it needs an exact treatment. Thus the degrees of freedom 
of the entire many body system are inexorably entangled 
with each other and therefore, it is no longer meaningful 
to separate what is a system from what is a bath. (5) 
Finally, a related point to (4) is in connection with the es- 
sential quantum nature of diamagnetism. As has been ar- 
gued by Jayannavar and Kumar |28j , not only is there no 
classical diamagnetism — due to the Bohr- Van Leeuwen 
theorem — there is no dissipative classical diamagnetism 
cither. Thus, the nonequilibrium, classical diamagne- 
tization relaxes to zero, a damping- independent result. 
The same is true for the classical damped harmonic os- 
cillator. In that case the time-dependent probability dis- 



tribution for the underlying Ornstcin-Uhlenbeck-process 
|29j relaxes to the equilibrium Boltzmann distribution, 
free of damping, even though the system-bath coupling 
is treated exactly through the classical Langevin equa- 
tions |3dj |. Therefore, we emphasize once again that the 
appearance of damping terms in equilibrium answers, as 
discussed under points (4) and (1), is an intrinsically non- 
classical aspect. 
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